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Abstract 

We study co-rotational wave maps from (3 + 1) -Minkowski space to the three-sphere S^. It is 
known that there exists a countable family {/n} of self-similar solutions. We investigate their stability 
under linear perturbations by operator theoretic methods. To this end we study the spectra of the 
perturbation operators, prove well-posedness of the corresponding linear Cauchy problem and deduce a 
growth estimate for solutions. Finally, we study perturbations of the limiting solution which is obtained 
from fn by letting n ^ oo. 

1 Introduction and Definition of the Model 

Many nonlinear evolution equations show a behaviour known as "blow up": Although the initial data 
are regular, the solution becomes singular in finite time. Important questions in this respect are: Under 
what circumstances does the blow up occur? What is the mechanism that leads to singularity formation? 
What is the blow up rate and the shape of the profile before the singularity forms? Unfortunately, many 
interesting examples from physics (such as Einstein's equations) are much too complicated to be accessible 
to presently available analytical techniques. Thus, one has to rely on numerical studies and simplified 
"toy models" which share some features with the original equation but are easier to handle. 

The wave maps system provides an example of such a toy model which has been studied in recent 
times. Consider (n + 1) -dimensional Minkowski space (R"''^^,rj) with metric {rjfj,„) — diag(— 1, 1, . . . , 1), 
an m-dimensional Riemannian manifold {M,g) with metric g and a mapping $ : R"+^ M. Choosing 
local coordinates on M we denote the components of $ and the metric g by and gAB, respectively, 
where capital latin indices run from 1 to m. The action functional 5* for wave maps is defined by 

where integration is understood with respect to the ordinary Lebesgue measure on 8"+^ and Einstein's 
summation convention is in force (greek indices run from to n). As usual we write := gfjr- The 
Euler-Lagrange equations associated to this action are given by 

□$-4 + ,,''-^r'*sc($)(a^<E>^)(c».$^) = (1) 

where bc ~ \g^^ {dBgcD + dcgBD — dogBc) are the Christoffel symbols associated to the metric g 
on M and □ := ri'^"df_idi, denotes the wave operator. Eq. ((T} is known as the wave maps equation and it 
is a system of semilinear wave equations. 

We are interested in the Cauchy problem for eq. ([l]), i.e. we prescribe initial data <I>'*|t^o = '^o and 
dt'^^\t=o = ^1 where dt :~ do and study the future development. There exists a well-established local 
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well-posedness theory for the wave maps system which states that for any initial data , G H" x 
H"-^ (s > ^) there exists a T > and functions $^ G C([0, T], with dt't'-^ G C{[0,T], H"-'^) that 
solve eq. ((T} in the sense of distributions and satisfy {^'^ , dt^^)\t=o = {^oj^i) where H" := H''{R") 
denotes the usual fractional Sobolev space. Moreover, the time T can be chosen to depend continuously on 
the initial data and the mapping , $^ ) ^ ($^, 9*$^) : H" x H"-^ C([0, T], i?") x C([0, T], 
is continuous (see [6] and references therein for details). There are also many recent results addressing 
global issues, see e.g. [16] for a review. Of most interest for our purposes is Tao's work ^TSl which deals 
with spherical targets, i.e. M is an (m — l)-sphere for m, n > 2. Under this assumption a global smooth 
solution of the wave maps system eq. ([1} exists whenever the initial data are smooth and small in the 
homogeneous Sobolev space H"^^ x H"^^'^. 

The present work deals with co-rotational wave maps from (3 + l)-Minkowski space to the three- 
sphere S^. Choosing the usual spherical coordinates (t, r, 6, ip) on Minkowski space and standard coordi- 
nates {ip^ Q, (f>) on with metric g — dijP +sin^ ip{dQ'^ +sin^ Qd(j}^), we restrict ourselves to mappings of 
the form {t,r,0,ip) [il){t,r),9,ip) which are called co-rotational. Under these assumptions the system 
eq. reduces to the single semilinear wave equation 

2 , sin(2i/;) 

V'tt - Vat - -Va- + = 0. (2) 

We remark that Tao's result applies to this problem but there are earlier papers by Sideris 13' and 
Shatah, Tahvildar-Zadeh [12] which address global questions for this particular model. 



2 Self— Similar Solutions and Singularity Formation 

Note that eq. ([2]) is invariant under dilations, i.e. if -i/) is a solution, so is t/ja defined by x}jx(t,r) := 
ip{t/ X,r / X) for a A > 0. There exists a conserved energy associated to eq. Q given by 

roc 

E^{t)= / (rV?(t,r)+rVr(i,r)+2sin'(V'(t,r)))dr 
Jo 

that scales as E^^ = which means that eq. © is energy supercritical. It is widely believed that 
for energy supercritical equations blow up is possible since it is energetically favourable for solutions to 
shrink which might result in the formation of a singularity. In order to obtain an explicit example of 
a blow up solution it is reasonable to look for self-similar solutions ip{t,r) = f{jr~t) where T > is a 
parameter (the blow up time). Plugging this ansatz into eq. ([2]) yields 

/" + -/'-^i^-0 (3) 
r P (1 — P ) 

where ' := ^ and p := jrr^- We require the regularity conditions /(O) = and /(I) = |-. Shatah [11] 
has shown existence of a smooth solution /o of eq. Q which has later been found in closed form [T7] and 
it is given by /o(p) = 2arctan(p). Bizoii [2] has extended this result by proving existence of a countable 
family {/„ : n = 0, 1, . . . } of analytic solutions of eq. ([Sjl satisfying /n(0) = and /n(l) ~ f where n 
counts the number of intersections of /„ with the line ^ on the interval [0, 1). Furthermore, the solutions 
/„ converge to the constant foo = -f for n — > oo pointwise on (0, 1]. These self-similar wave maps provide 
explicit examples of solutions with smooth finite energy initial data (use appropriate cut-off functions 
and finite speed of propagation) that become singular when t T— . 

An obviously important question is whether blow up occurs for generic initial data or the self-similar 
examples are "exceptional cases". This issue has been addressed by Bizoh et. al. [^ using numerical 
techniques. Their results strongly suggest that the self-similar solution fo provides a universal blow up 
profile. More precise, they conjecture that there exists a large set of initial data that lead to blow up 
and the singularity formation takes place via fo, i.e. the time evolution approaches fo near the center 
r = as t ^ T-. 

We remark that this self-similar blow up behaviour is very different from singularity formation in 
the energy critical case of co-rotational wave maps from (2-|-l)-Minkowski space to the two-sphere. 
For this model it has been known for quite some time (see [TJ]) that the blow up rate cannot be self- 
similar. However, rigorous existence of blow up solutions has only recently been shown by Rodnianski 
and Sterbenz TOl as well as Krieger, Schlag, Tataru [?]■ 
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3 The Linear Stability Problem 



As mentioned, the self-similar solution /o seems to play an important role in dynamical time evolution. 
Thus, it has to be stable in a certain sense since otherwise it could not be found by numerical time 
evolutions starting from generic initial data. In order to study linear stability of the wave map /„ we 
introduce adapted coordinates (a, p) given by cr = — log \/ {T ~ t)^ — and p = jt—^ (cf. The 
introduction of similarity coordinates in the study of self-similar blow up is very common (see e.g. [8]). 
We refer to {a, p) as hyperbolic coordinates since the lines a — const are hyperbolae in a spacetime 
diagram. Note that the coordinate system {a, p) is only defined for r < T — t and thus, it covers the 
interior of the backward lightcone of the blow up point {t, r) — (T, 0). Furthermore, (a, p) is an orthogonal 
coordinate system, i.e. the vectorfields dc and dp are orthogonal with respect to the Minkowski metric. 
Moreover, the blow up time t = T is shifted to infinity, i.e. t — > T— corresponds to cr ^ oo. The valid 
range is ct £ R and p G [0, 1). Eq. ((2)1 transforms into 

2x2, 2(l-p2)2 , (1 - p2) sin(2?/;) ^ 

p p2 

Plugging the ansatz 'ip[a,p) = /n(p) + <l>{i7,p) into eq. Q and linearizing around /„ we arrive at the 
evolution equation 

^ 9^ n 2^2^ 2(1 -p^)^ , , 2(1 - p^) cos(2/„) ^ 

(paa — 2(pa — (1 — p ) (ppp (f>p H 5 (p = (5) 

p p^ 

for the perturbation which governs the linearized flow around the wave map /„. The wave map /„ 
is considered to be linearly stable if there do not exist solutions of eq. ((SJ that grow (with respect to 
a suitable norm) as a increases. We will make this more precise below. However, a direct computation 
shows that (j)^ defined by (l)n{o,p) := e'^ p^l — p^fn{p) solves eq. ([5]). This instability is related to the 
time translation symmetry of the original problem. The reason for this is the fact that /„ is not a single 
solution but a one-parameter family of solutions depending on the blow up time T. We will refer to this 
instability as the gauge instability. Thus, without further considerations the best we can expect is to 
prove a result that excludes the existence of solutions of eq. (O with /„ = /o that grow faster than this 
gauge instability. In what follows we will address this issue by using an operator theoretic approach. 



4 Operator Formulation 

In order to simplify the problem further we define a new unknown by 0((j, p) := e~''(j){a,p) to remove 
the first-order term —2(j)cr- Eq. ([Sjl transforms into 

- (1 - pY0.. - + - ^"^ -^^ = 0. (6) 



4.1 The Operator An 

We define the formal Sturm-Liouville differential expression a„ on the interval (0, 1) by 

a„u ~ —(—(pu'y + q„u) 
w 

where w{p) :— -^j^-^^, p{p) ■= p^ and gn(p) )^°^'^^Mp))-p ^ Thus, a„ represents the spatial 

differential operator in eq. ([6)|. We define the weighted Hilbert space H :— 1/55,(0,1) where L55,(0, 1) 
consists of (equivalence classes of) functions which are square integrable with respect to the weight w, 
i.e. the inner product {■\-)h on H is given by 

{u\v)h = / u{p)v{p)w{p)dp 
Jo 

and as usual we write || ■ ||h ~ \/ {■\-)h- We set 2?(An) :— {u £ H : u,pu' G ylCioc(0, l),anU G H} and 
AnU :— anU for u £ T>{An). 



3 



Lemma 4.1. The operator An : T>{An) C H ^ H is self-adjoint. 

Proof. Consider solutions of a„u — 0. Since the coefficients of tliis linear ordinary differential equation 
are meromorphic and p = as well as p = 1 are regular singular points, we can apply Frobenius' method 
to obtain series representations for solutions around p = and p = 1. Around p — the Frobenius 
indices are —2 and 1, i.e. only one nontrivial solution belongs to H . According to the Weyl alternative 
we conclude that p = is in the limit-point case (cf. e.g. [20]). Similarly, around p = 1 the indices are 
both ^ and thus, neither of the two linearly independent solutions belongs to H which shows that p = 1 
is in the limit-point case as well, ft follows that the maximal operator An associated to a„ is self-adjoint 

(cf. ED]). □ 
4.2 Factorization of Aq 

Very similar to the energy critical case [10] there exists a factorization of the linearized operator that 
immediately implies its nonnegativity. To see this, note first that the potential term in eq. ([5} for the 
ground state fo can be rewritten according to 

cos(2/o(p)) = cos(4arctanp) = ^ j" Z' . 

(1 + p'')'' 

Now we define two first-order operators B and 13 on H hy T>{B) := T>{13) := C^(0, 1) as well as 

Bu{p) ~ (1 - p^)u'(p) - -7T-^75tw(p) 



and 



g _ 2 

Bu{p) := -(1 - p'^)u'{p) 77— -2T"W- 



A simple integration by parts shows that {Bu\v)h = {u\Bv)h for all u,v £ C^(0, 1) and thus, B and 
B are adjoint to each other. Furthermore, by a direct computation one verifies that Aqu = BBu for 
all u G C^(0, 1) which is the desired factorization. We immediately obtain (Aou\u)h ~ {BBu\u)h = 
{Bu\Bu)h > for all u e C^(0, 1). Since Ao is limit-point at both endpoints it follows that C^(0, 1) 
is a core for Aq and the estimate {Aou\u) > extends to all u £ 'D(Ao) by a density argument. Thus, 
Ao is nonnegative and we have a{Ao) C [0, oo). 

4.3 The Spectrum of 

Proposition 4.1. The point spectrum ap{An) of An consists of exactly n negative eigenvalues, i.e. 
(jp{An) — {Ai,...,A„ : > Ai > ■■■ > A„}. Furthermore, the essential spectrum ae{An) is given by 
CT^(An) = [0,oo). 

Proof. Consider the equation (A — an)u = 0. The Frobenius indices around p = 1 are ii^^. Thus, for 
A > neither of the two linearly independent solutions belongs to H which shows (Tp(yl„) C (— oo,0). 
Furthermore, we conclude ae{An) = [0, oo) by |18] . p. 165, Theorem 11.5c. 

As already discussed, the function On defined by 6n{p) := P\/l — P^fLip) is a solution of a„M — 0. 
Bizon [2] has shown that fn{p) £ (0,7r) for p £ (0, 1] and thus, if /„ > .|, it can only have one maximum 
since 

„ sin(2/„(p)) 
p2(l - p2) 

for fn{p) ~ (cf. eq. (O). Similarly, if /„ < |-, it can only have one minimum. Hence, the number of 
extrema of /„ equals the number of intersections of /„ with the line |- on (0, 1) and therefore, dn has 
exactly n zeros on (0, 1). 

Applying the transformation v{x) :— sinh(x)it(tanha;), the equation (A — a„)u — reduces to the 
eigenvalue problem 

Xv{x) + v"{x) - ^v{x) ~ Qn{x)v{x) = (7) 

X'' 

where 

^ 2cos(2/„(tanha:)) _ 2_ 

sinh^ X x^ 
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Consequently, the function ifin{x) := sinh(a;)S„(tanh2;) is a solution of eq. ((T} with A = which has 
exactly n zeros on a; £ (0, oo). Observe that the mapping u i— > i; is an isometry from 1/^(0, 1) to L^{0, oo) 
and thus, the study of An is equivalent to the spectral problem for the £ = 1 radial Schrodinger operator 
on L^(0, oo) with potential Qn- Note that Q„ £ I/°°(0, oo) and from the Frobenius analysis of a„ it 
follows that the asymptotic behaviour of nontrivial solutions of eq. ([7} is the same as if Q„ G (0, oo) . 
Thus, the results of [9J, XIII. 3. B apply and we conclude that A„ has exactly n eigenvalues below the 
essential spectrum. □ 

Remark 4.1. This observation has already been made by Bizoh T. Note further that 9n is not an 
eigenfunction since it does not belong to H\ 

Remark 4.2. Negative eigenvalues of An lead to unstable mode solutions of eq. ([5]), i.e. <j){a,p) :— 
g(i+v'^)<T^j-^'j jg ^ solution of eq. ([5]) if A < belongs to ap{An) and v is the associated eigenfunction. 
Thus, the n-th self-similar solution has exactly n unstable modes. 



5 Well— Posedness of the Linearized Problem and Growth 
of Solutions 

We consider the linear Cauchy problem 

u{a) + Aou{a) = for cr > , . 

u(0) = uo,u{0) = til ^ ' 

for a function u : [0, oo) H and uo,ui € 'D{Ao) where ' :— -j^ and differentiation with respect to a is 
understood in the strong sense. Eq. ^ is an operator formulation of eq. ([SJ with /„ = /o and thus, it 
describes the (rescaled) linear flow around the self-similar wave map fo- 

According to the spectral theorem for self-adjoint operators (see e.g. [l^) there exists a spectral 
family {E{X) : A £ K} of orthogonal projections associated to Aq such that for u £ T>{Ao) and v G H we 
have 



{Aou\v)h = J idud^u,v 

where Hu.v ■— j{n /Xu is the unique measure defined on the Borel- 

(T-algebra on R that satisfies /iu((Ai,A2]) = ||(i?(A2) — E{Xi))u\\^ for u £ H. Moreover, with the help 
of the spectral family {E{X)} it is possible to define functions of self-adjoint operators. Let / : R — > C 
be measureable and set V{f(Ao)) :— {u £ H : J\f\^dfiu < oo}, {f{Ao)u\v)H :~ J fd^u,v for any 
u £ 'D{f{Ao)) and v £ H. Note that in fact only the values of / on the spectrum of Aq contribute. 
This yields a well-defined linear operator f{Ao) : T){f{Ao)) C H ^ H which behaves nicely (see [19] for 
details and a functional calculus). 

Proposition 5.1. Let uo,ui £ I?(j4o). Then, the Cauchy problem eq. (8^ has a unique solution u : 
[0, oo) — > V{Ao) C H and the real-valued function a i-^ \\Af/'^u{a)\\^ + \\u{a)\\'jj is constant for all 
o- > 0. 

Proof. We set 

, , _ J sin(crA)A"^ for A > 
ga[X) I ^fo, ;,<o 

l/2 1/2 1/2 

and define u{a) := cos((tAq )uo + gaiA^ )ui via the functional calculus. Note that Aq is well- 
defined since a{Ao) — [0,oo). Applying Lebesgue's theorem on dominated convergence one immediately 
concludes that ti((T) + Aou{a) = and clearly, u satisfies the initial conditions. Furthermore, it follows 
easily that u{a) and 'u,{a) stay in ©(Ao) for all a > 0. A direct computation (functional calculus) yields 
\\Al^\{a)\\j, + \\u{a)\\ \ = \\A\ ^moIIh + ||iii||l/ which implies uniqueness as well. □ 

Note that (uo,tti) ||ylQ''^ito|||/ + ||iii||lf defines a norm on T){A^q'^) x H since ^ ap{A^J'^) and 
hence, there do not exist growing solutions of eq. ^ which is the desired result. Thus, in view of the 
transformation (j) i-^ (f> we have shown that there do not exist linear perturbations of fo that grow faster 
than the gauge instability. 
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6 Discussion 



We give a heuristic explanation why one cannot go any further in this formulation. This is most easily 
seen within a first-order formulation. Let u : [0, oo) — > _ff be the unique solution of eq. ^ with initial 
data uo,ui € 2?(Ao) and set </'(o", ■) := e''u{a). Then, (f) satisfies the first-order equation 

where u{a) := (<^(cr, •), 4><t{cf, •)) and 




4> (resp. u) describes the linear flow around the self-similar wave map /o. Consider the operator 
L on the Banach space X :— 'D{Ag'''^) x H with norm || (ui , 1*2)1] ^ := ||ui||^i/2 + ll'falllf and ©(L) := 

'D{Ao) X 'D{Aq^^) where |j • || 1/2 denotes the graph norm of A^^. Then, it is easily seen that the spectrum 

of L can be calculated from the spectrum of Aq and it is given by o"(L) — {z € C : —{z — 1)^ £ a{Ao)} 
which shows that o-(L) — {z £ C : Ke{z) — 1}. We have already discussed the gauge instability which 
corresponds to the point 1 £ f^iL). Thus, the origin of the point 1 G (^{L) is well-understood and in 
fact, we are only interested in stability of /o "modulo this gauge instability". If 1 G o-{L) would be an 
isolated point one could define a projection operator to remove this instability. However, the rest of the 
spectrum makes this impossible and in what follows we give an intuitive explanation where it emanates 
from. The reason for this is a serious defect of the coordinate system {o,p). As already mentioned, it 
covers only the interior of the backward lightcone of the point (T, 0). Thus, an outgoing wave packet 
can never leave the backward lightcone but cumulates near p = 1 as cr increases. This corresponds to 
an apparent instability which shows up in the spectrum of L. Unfortunately, one cannot go around this 
difficulty unless one changes coordinates which destroys the self-adjoint character of the problem. Thus, 
the result that there do not exist linear perturbations of /o that grow faster than the gauge instability is 
the best one can have in this formulation but still, it strongly supports the conjecture of stability of the 
wave map /o. 



7 The Operators An for n ^ oo 

The self-similar wave maps /„ can be constructed numerically using a shooting and matching technique 
[1]. By the same method one may calculate the point spectrum of An (cf. [4]). Results for Ai, . . . , A4, 
are given in Table [T] where fij := \/— Aj. 

7.1 The Operator 

We emphasize the remarkable fact that the horizontal rows in Table [1] seem to converge. The question is 
if one can construct an operator Aoo such that the point spectrum of Aoo is the limit of the point spectra 
of An for n 00. Since /„ /oo = -f pointwise on (0, 1], the obvious strategy is to investigate linear 
perturbations around foo described by eq. ((6} with /„ = /oo. Thus, we consider the Sturm-Liouville 
expression 

flooW := —{—{pu'Y + qoou) 
w 

where ?oo(p) := f^^'^iyi ■ Again, we use Frobenius' method to classify the endpoints. The behaviour of 
solutions of flooW = around p = 1 is the same as before for a„ and thus, p = 1 is in the limit-point 
case again. However, around p = the Frobenius indices are and thus, all solutions of OooW = 

belong to H near p — 0. According to the Weyl alternative the point p = is in the limit-circle case 
and hence, we have to specify a boundary condition of the form [u, x]p(0) = for a fixed function x 
in order to obtain a self-adjoint operator (cf. [20)1. The question is how to choose x- Bizoh has 
studied an analogous problem for a Yang-Mills system. Following his proposal we choose x in such a 
way that UooU ~ has a nontrivial solution just like the gauge instability for a„. Thus, we denote by 
X a real-valued nontrivial function satisfying aooX = with asymptotic behaviour x(p) ~ (1 — p)^^'^ for 
p — > 1. Such a function exists by Frobenius' method and it is unique up to constant multiples since the 
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second linearly independent solution around p = 1 contains a log term. Using an appropriate smooth 
cut -off function we construct a x that coincides with x on (0, i) and satisfies Xj'JooX G -ff as well as 
X,PX' G ACioc(0, 1). Then, we define 

f (^oo) --{ue H : u,pu G ACioc(0, l),acx,tt £ [m,x1p(0) = 0} 

and AaoU '■= a-ooU for u G ©(Aoo). Invoking Sturm-Liouville theory ^0] we conclude that Aoo ■ T){Aaa) C 
H ^ H IS self-adjoint. 

7.2 The Point Spectrum of 

Since the asymptotic behaviour of solutions of (A— Ooo)'" = around p = 1 is the same as for (A— a„)it = 0, 
we conclude that ap{Aoo) C (— oo,0). Thus, we restrict ourselves to A < 0. The eigenvalue equation 
(A — aoo)u = is a linear ordinary differential equation of second order with four regular singular points 
p = 0, ±1,00. However, applying the transformation p z :— fp' we can reduce this number by 1 and 
we are left with the three regular singular points z = 0,1, oo. Thus, solutions are given in terms of 
hypergeometric functions (cf. t5j). We define a new unknown 

v[z) := z^il — z)^u[^/z) 

where a := —^(—1 + i\/7) and /3 ~ ^1(1 + V^-^)- The equation (A — aoo)u = transforms into 

z{l — z)v" + [c — (a + 6 + l)z\v' — abv — 

where a := \{1 + 2^/^ + i^/7), b := 1(3 + 2^^ + ^^) and c := 1 + The solution t;i(-. A) defined 
by Vi{z, A) := 2-Fi(a, b; a + b + 1 — c; 1 — z) corresponds to the solution of (A — aoo)w = which belongs 
to H near p = 1. vi can be written as a linear combination of the form 

Vi {z, A) oc 2-Fi (a, b; c; z) + m{X)z^^'^2Fi (a + 1 — c, 6 + 1 — c; 2 — c; 2) 

where the connection coefficient m(A) can be given in terms of the F-function [S] and reads 

r(a + 1 - c)r(6 + 1 - c)r(c - 1) 



m(A) = 



r(a)r(6)r(i 



Thus, the boundary condition at p = translates into m(A) — m(0). This transcendental equation can 
be solved numerically and the results are given in Table [T] 







A2 


^3 




^00 


Ml 


5.333625 


5.304 


5.30 


5.3 


5.3009 


^J■2 




58.0701 


57.68 


57.6 


57.637 








625 


620 


619.61 



Table 1: Approximate eigenvalues 
Thus, as conjectured, the point spectrum of Aoo seems to be the limit of the point spectra of An 



8 Generalization to higher equivariance indices 

Finally, we consider more general equivariant wave maps described by the equation 



2 , £(£+1) sm(2ih) 
r 2 r-' 



(9) 



where G N is the equivariance index (see e.g. [12] for a derivation of this equation). Self-similar 
solutions •il){t,r) — /(jt^) of eq. (|9} satisfy 



_;^+l) sin(2/) 



p2(l 
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and again, we require /(O) — and /(I) = ^ as regularity conditions. Bizori's proof of the existence of 
a countable set {fn,e : n = 0, 1, 2 . . . } of self-similar solutions carries over to this more general situation 
(see In fact, the solutions f^^i have analogous properties for all i. The index n counts the number 
of intersections of f„,e with the line ^ in the interval p £ [0, 1) and we have fn,e{p) G (0, tt) for all n and 
p G (0, 1]. Furthermore, the asymptotic behaviour of for p — > is fn,e{p) ~ P^- 
Linear perturbations (f> around f„,e satisfy the evolution equation 

rh OrK n ^ii^Zl!^ , ^(^ + 1) ( 1 - p" ) COs(2/„., ) ^ 

4>an — 20CT — (1 — p ) (ppp 4>P H 5 9 = 

p p^ 

and analogous to the case ^ = 1 we define the Sturm-Liouville expressions a^^i by 

a„,£U := —(—{pu')' + qn,tu) 
w 

where w{p) ■- jr^iyi, p(p) — and 

^(^+l)(l-p2)cos(2/„,,(p))-p2 



qnAp) ■= 



(1 



The Frobenius indices for the equation a„,£it = are — 1,^} at p = and {i, i} at p = 1 and 
thus, both endpoints are in the limit-point case. It follows that the associated maximal operator An,e 
on H = 1/^(0, 1) is self-adjoint. Again, the study of A„^e is equivalent to the spectral problem for the 
radial Schrodinger operator and thus, the rest of the analysis remains unchanged. Hence, we have the 
following theorem. 

Theorem 8.1. The spectrum of A„^e consists of a continuous part ae{An.t) — [0,00) and exactly n 
negative eigenvalues. As a consequence, Aq^( is nonnegative, the Cauchy problem 

ilia) + Ao,eu{a) = for a > 
m(0) = uo,u{0) = iti 

foruo,u\ G T>{An^i) is well-posed and the real-valued function a 1-^ \\AQ^fu{a)\\^ + \\u{a)\\'jj is constant 
for all a > 0. 
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